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In this paper, we present new progress on the study of the symmetric extension criterion for 
separability. First, we show that a perturbation of order 0{1/N) is sufficient and, in general, 
necessary to destroy the entanglement of any state admitting an A'' Bose symmetric extension. On 
the other hand, the minimum amount of local noise necessary to induce separability on states arising 
from Bose symmetric extensions with Positive Partial Transpose (PPT) decreases at least as fast 
as 0{\/N^). From these results, we derive upper bounds on the time and space complexity of 
the weak membership problem of separability when attacked via algorithms that search for PPT 
symmetric extensions. Finally, we show how to estimate the error we incur when we approximate 
the set of separable states by the set of (PPT) A^'-extendable quantum states in order to compute the 
maximum average fidelity in pure state estimation problems, the maximal output purity of quantum 
channels, and the geometric measure of entanglement. 

PACS numbers: 



I. INTRODUCTION 

The separability problem, that is, the problem to de- 
termine whether a given quantum state is separable or 
entangled, is one of the most fundamental problems in 
Entanglement Theory [l[ . Starting from the famous PPT 
(Positive Partial Transpose) criterion [5], nowadays we 
have an enormous number of different separability crite- 
ria to choose from (see the citation lists of review papers 
in this topic [J, i, i, S i, 0, H)- Among all known 
separability criteria, those based on "symmetric exten- 
sions" and "PPT symmetric extensions" (i.e., symmetric 
extensions with an additional PPT constraint), as con- 
ceived by Doherty et al. S are considered to be 
among the most powerful fa]. These criteria rely on the 
fact that any set of A^-symmetrically extendable states 
(PPT or not) converges to a set of separable states in 
the limit of TV oo, as first noticed by Raggio and 
Werner [Til , [l^ , although it also follows from the Quan- 
tum de Finetti theorem Since both the set of A^- 
symmetrically extendable states and the set of A^-PPT 
symmetrically extendable states can be characterized by 
Semidefinite Programming p^ . a well-known optimiza- 
tion problem for which many free solvers are available 
(like the MATLAB toolbox SeDuMi[ll]), these tests are 
not only powerful, but also easy to implement. This ex- 
plains why, over all known numerical methods, the algo- 
rithms created by Doherty, Parrilo and Spedalieri (DPS) 
are the most popular in the Quantum Information com- 
munity (notice, however, that there exist other methods 
for entanglement detection based on Semidefinite Pro- 
gramming besides the DPS criterion [l^[l3|)- 

This family of schemes has, though, an important 
drawback: in this approach, in order to conclude that 
a given state p is entangled, it is enough to find an A^ 
such that p does not belong to the set of A^-(PPT) sym- 
metric extendable states. On the other hand, in order 
to show that a given state is separable, we would have 



to prove that it admits an A^-(PPT) symmetric exten- 
sion for all natural numbers N . The DPS method then 
becomes useless: since we always operate under finite 
time and memory constraints, all we can do in practice 
is to check for the existence of A^-(PPT) symmetric ex- 
tensions for A^ less or equal than some finite number A'q. 
If the state p under analysis happened to admit an A^o 
(PPT) symmetric extension, we could thus not conclude 
anything about its separability. 

Hulke and Bruss [l^ tried to solve the issue by pro- 
viding a complementary criterion designed to detect sep- 
arability instead of entanglement, to be implemented at 
the same time as the DPS criterion. Unfortunately, the 
time complexity of that other method scales superexpo- 
nentially with the dimension of the subsystems involved 
0. The reduced speed of convergence of the resulting 
two-way algorithm (much smaller than that of the DPS 
criterion) thus makes it unsuitable to study quantum cor- 
relations in high dimensional systems. 

Besides, there is a more elegant way to approach the 
problem. 

In a recent work, loannou observed that, even if a state 
happens to have an A^o-(PPT) symmetric extension, we 
can at least bound the distance between such state and 
the set of separable states in terms of A^o • In the lan- 
guage of Computer Science, this means that the "trun- 
cated" DPS criterion allows to solve an instance of an 
approximate separability problem, the weak membership 
problem of separability {WMEM{S)). loannou therefore 
provided an upper bound on the full time-complexity of 
the algorithm for WMEM based on symmetric extension 
criteria. 

But even after loannou's work, an open question re- 
mains to be solved. The PPT symmetric extension cri- 
terion is considered to be stronger than the symmetric 
extension criterion [lo| . By definition, it is actually 
at least as strong as the symmetric extension criterion in 
the sense that a A^-PPT symmetrically extendable state 
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is TV-symmetrically extendable. However, so far, there 
are no results that quantify how strong the additional 
PPT constraint makes the DPS criterion. In particular, 
since the additional PPT constraint increases quadrati- 
cally the size of the matrices that define the Semidefinite 
Programming problem, there still remains the possibil- 
ity that the PPT criterion just makes the DPS algorithm 
slower for WMEM{S). In order to make this point clear, a 
similar analysis as loannou's should be done for the PPT- 
symmetric extension criteria. Since loannou's analysis is 
based on the finite quantum de Finetti theorem H^l 
and there exists no similar theorem for states satisfying 
the PPT constraint, there is no straightforward exten- 
sion of loannou's work to the PPT symmetric extension 
criterion. 

In this paper, by analyzing these criteria in more de- 
tail, we extend loannou's result to account for the PPT 
condition. 

The structure of this article is as follows: in Section 
im we will give the reader a detailed explanation of the 
DPS criterion and introduce the basic notation that will 
be used in the paper. Then we will move on to present 
the main result of this article, namely, an upper bound 
on the amount of noise needed to make the DPS states 
separable. This will allow us to compute upper bounds 
on the entanglement robustness of these states, and on 
their distance to the set of separable states. We will 
also briefly discuss how close our bounds are to being 
optimal. In Section IIV| we will use the previous results 
to analyze the computational complexity of solving the 
weak membership problem of separability through the 
DPS criterion. In particular, we will show that the PPT 
constraint in the DPS criterion reduces the dominant 
factor of the upper bound on the time complexity from 
[ki/Sf"^'^ to (k2/ 5)'^'^'^ , where 5 is the accuracy parame- 
ter of WMEM{S). In Section |V] we wiU bound the speed 
of convergence of the DPS criterion when applied to com- 
pute the optimal fidelity in state estimation problems, 
the output purity of quantum channels and the geometric 
entanglement of arbitrary states. There we will perform 
some numerical tests to have a grasp at the actual speed 
of convergence of the DPS criterion, as opposed to our 
analytical upper bounds on it. In Sections IVIl IVIII we 
will give the proof of the main theorem and explain how 
it can be extended to deal with the multipartite case. 
Afterwards, we will also show a very simple method to 
bound the entanglement of general PPT states. Finally, 
Section [IXI will present our conclusions. 



II. THE DPS CRITERION 

The Doherty-Parrilo-Spedalieri (DPS) criterion for en- 
tanglement detection |^ is a numerical algorithm that, 
combining the aforementioned results [111, [3, [l^ on N- 
extendibility with convex optimization methods, allows 
to characterize the set S of separable operators up to ar- 
bitrary precision. The criterion arises from the following 



observation: if Kab G then, by definition, it belongs 
to the cone of bipartite product states, i.e., 

Aab = ^Pi\ui){u.^ (g) \v^){v^\, (1) 

i 

with Pi > for all i. 

Once this decomposition is known, we can define a 
uniparametric family of operators A^^jv g B{J-La®'H%'^) 
by tensoring N times the last part: 

Kab- = Y.P^\u,){u,\ ® \v,){v,f>'' . (2) 

i 

Let us study the properties of the newly defined opera- 
tors: first of all, from the above definition it is clear that 
they are all positive semidefinite. Also, from ([2]) it can be 
seen that tracing out the last iV— 1 systems we recover the 
initial operator, i.e., tiB'^-^i^AB") = ^abj and that the 
last N systems are invariant under the action of the per- 
mutation group. Finally, when viewed as an A^+l-partite 
system, A^^w is multiseparable, and therefore must re- 
main positive semidefinite under the partial transposition 
of any bipartition of these systems. 

For simplicity, we will incorporate all these properties 
in a single definition: 

Definition 1. Bose symmetric extensions (BSE) 
Let h-AB G B{TLa ®'Hb) be a non-negative operator. We 
will say that Aab- € ^O^A ® Ti,^^) is an N Bose sym- 
metric extension (BSE) of Kab iff- 

1. Kabn > 0. 

2. trBN-i{KAB-) = Aab- 

3. A^^N is Bose symmetric, i.e., KAB-{^A®Psym) = 
A^^N, where P^ym denotes the symmetric projec- 
tor of N particles. 

In case Kab- is PPT with respect to all or some of its 
bipartitions AB^ \B^~^ , we will call it a PPT Bose sym- 
metric extension (PPT BSE) of Ayi^. 

From what we have seen, it is clear that, if Kab is a 
separable operator, then there exists an N (PPT) BSE of 
Kab for any TV. Since (PPT) Bose symmetric extensions 
are defined through linear matrix inequalities, the prob- 
lem of determining whether a given state Kab admits 
one or not can be cast as a semidefinite program (SDP) 
and therefore can be solved efficiently for fixed N 
and varying dimensions. The DPS criterion consists pre- 
cisely in, given an operator A^^ whose separability is at 
stake, check for the existence of N (PPT) Bose symmet- 
ric extensions for different values of N. 

A hierarchy of separability tests arises then naturally: 
if some operator Kab does not admit a (PPT) Bose sym- 
metric extension for some N (i.e., it does not pass the 
iV*'* test), then it has to be entangled. If, on the con- 
trary, such extension exists, then we would go for the 
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(TV + 1)*'' test, that is, we would search for A?^ + 1 (PPT) 
Bose symmetric extensions of A^^ . This last test would 
be in general more restrictive than the previous one, since 
for any N + 1 (PPT) Bose symmetric extension A^^n+i 
of Aab we can obtain an N (PPT) Bose symmetric ex- 
tension by tracin g ou t the last system. 

Doherty et al. [13] showed that the previous hierarchy 
completely characterizes the set of separable operators, 
in the sense that for any entangled positive operator Aab 
there exists an N such that Aab does not admit an N 
Bose symmetric extension. 

We will now introduce a notation that will be used for 
the rest of the article: will denote the cone of all 
bipartite operators that have an N BSE, and Sp will 
be understood as the set of all unnormalized quantum 
states that admit an N BSE that is PPT with respect to 
the bipartition AB^^^^'^ ^j^l^/^l , case we also demand 
normalization, we will be dealing with the sets of states 
S^,Sp. The elements of the previous four sets will be 
called A^-(PPT) symmetrically extendable operators, or 
states, if normalized, or just DPS operators or states. 
Our previous discussion can then be summarized as 



S^D S^D S^D ... D S, 
Sp ^ Sp ^ Sp ^ ... ZD »5*, 



lim 



N- 



.s^,s!^ 



s. 



(3) 



Note that S^ = S^ (5"^ = 5,^) is the set of all positive 
semidefinite operators (states). 

Before ending this section, we would like to point out 
one additional fact. As we already explained in the in- 
troduction, when we use the DPS criterion in practice, 
it is not possible to conclude with certainty that a given 
state is separable. However, in the PPT case, by check- 
ing some rank constraints on the density matrices output 
by the computer, we can sometimes conclude separabil- 
ity in a finite number of steps. In that case, we will say 
that the PPT BSE presents a rank loop. We will make 
use of rank loops in Section |V] in order to estimate the 
accuracy of our upper bounds on the error we introduce 
when we perform linear optimizations over the sets S^ 
or Sp instead of S in state estimation problems. A de- 
tailed explanation of this criterion for optimality can be 
found in Appendix IB] 



III. CHARACTERIZATION OF 5"^ AND Sp 

We have seen that the sequences of sets (S^), (Sp) 
tend to the set S in the limit N ^ oo. Intuitively, this 
means that, for iV >> 1, any state pab belonging to one 
of these sets must be either separable, or, at least, very 
close to a separable state. It seems thus plausible that 
the little entanglement such states may possess could be 
destroyed by some very attenuated local noise. One of 
the most simple noise models one can think of is depo- 
larization, where a quantum state is turned into white 



noise with probability p. The action of the depolarizing 
channel fi'^^-' over some state p € B{Ti) is given by 



n(f)(p) = (i-p)p + p-. 



(4) 



where d is the dimension of the Hilbert space Ti.. Given 
any bipartite quantum state pab, shared by Alice and 
Bob, we could thus define its critical disentangling prob- 
ability Pc{pab) as the minimum probability with which 
one of the parties, say Bob, would have to prepare the 
maximally mixed state in his subsystem in order to dis- 
entangle it from Alice's. That is, 



Pc{pab) = min{p : 1^ 



n^^\pAB)eS}. (5) 



Similarly, we can define the critical disentangling proba- 
bility of a set of states W as the maximum of all Pc{p) 
for all p S W. Clearly, Pc < ^ for all states, although this 
bound can be greatly improved if the dimensionality of 
Bob's system is small, as we shall see. 

In this section, we will give upper bounds on this crit- 
ical probability valid for any state in S'^ (or Sp). Then, 
by means of these results, we will provide several upper 
bounds on the speed of convergence of S'^ and Sp to S. 

Before proceeding, though, a remark on notation: in 
this article, we will be mainly concerned with linear op- 
erators or quantum states acting over a bipartite Hilbert 
space Ha (X" Hb , and all the formulas and bounds that 
we will derive in this section and the following three will 
involve the dimension of the Hilbert space Ti.B where 
the symmetric extensions are to be made. For the 
sake of clarity, we will therefore introduce the notation 

d dim Hb ■ 

The following theorems will play a key role in deriving 
most of the results of this paper. 



Theorem 2. 



Pc(^^) < 



N + d 



(6) 



In other words: for any operator Aab G S , the positive 
semidefinite operator 



A 



N 



AB 



-A 



1 



N + d"^^ 'N + d' 



-A 



is separable. 



Theorem 3. Define g^ (or g^^j^^ in case d is ambiguous) 
as 

gN = min{l — x : ^jv/2+i^(^) ~ ^} f'^''^ ^ even, 

min{l - X : P^N'+^y^i^) = 0} for N odd, (8) 

with Pji"'^^ (x) being the Jacobi Polynomials fsJi l. 



A'iSl llB 



(7) 
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Then, 



Corollary 4. Any p £ S satisfies 



9N- 



(9) 



That is, for any Aab G Sp , the positive semidefinite 
operator 



Aab = (1- 



is separable. 



2{d-l) 



on) Aab - 



2{d~l) 



qnAa^Ib (10) 



The proof of these two theorems is given in Section IVIl 
where a separable decomposition for the states ([7]), (fTO]) 
is also provided. Also, it is worth mentioning that, in 
both cases, A^i = Aa- 

Notice that, in Theorem [31 (Jat is defined in terms of 
the greatest root of Jacobi polynomials. The properties 
of the roots of Jacobi polynomials have been studied for 
quite time plj . This allows us to derive an expression 
for the asymptotic behavior oi g^: 



9n 



Jd-2,1 

N 



for N » 1, 



(i+ 1.856^1/3 + 0(^-1/3) 



N 

for iV > d > 1, 



(11) 



(12) 



where jn,i is the first positive zero of the Bessel function 
Jn{y)- 

How far can then the states in , be from the 
set S of separable states? A way to answer this question 
could be to bound the maximum possible entanglement 
of such states. 

The robustness of entanglement of a state p is defined 
as the minimum amount of separable noise needed to 
destroy the entanglement of such a state [121 '■ 

R{p) =^ mm{A : 3a G 5, s.t. e S}. (13) 

The robustness of entanglement is also an upper bound 
on the global robustness of entanglement Rg{p) [13]' 
fined by allowing a to be an arbitrary normalized quan- 
tum state in the above expression. And the global ro- 
bustness of entanglement is, in turn, lower bounded by 
several other entanglement measures, like the negativity, 
the geometric measure of entanglement and the relative 
entropy of entanglement [H, [H, [13, [H, 26]. Any non 
trivial upper bound on the entanglement robustness of 
the states in and Sp could thus retrieve a lot of in- 
formation. 

The following corollary follows straightforwardly from 
theorems [5] and [31 



R{p)< 



d-l 



N 



Similarly, any p £ Sp satisfies 

9n 



R{p)< 



(14) 



(15) 



To see why, suppose that p is normalized and use for- 
mulas (O, pUjl to express p (i.e., A^^) in each case 
as a convex sum of the non negative operators p and 
a = 3zj(pA ®^B — p)- Then, notice that, since pA = Pa 
and p is separable, then a must also be a separable op- 
erator [6lj 

Theorems [H and [31 also allow to obtain bounds on the 
distance between the states in pab & S^ , Sp and the set 
of separable states S. 

Corollary 5. For any p S S^ , there exist p E S such 
that 



\P-P\\i < 

]P- P]|oo < 

\p- pWf = 



2{d-l) 
N + d~l' 

d-l 
N + d-l 
d 



(16) 
(17) 



N + d 



Mp^)-^, (18) 



where \\ ■ \\ ■ ]|oo o,nd \\ ■ \\p are the trace, the operator 
and the Frobenius norm, respectively. 

Similarly, for any p g Bp ( and N > 2), there exists a 
state p G S such that 



\p- P\\i < 9N, 
]p-/5]|oo < 9n/2, 

I ~|| d9N 
\P- P\\ 



2d -2 



tr{p^) 



HPa) 



(19) 
(20) 

(21) 



Proof. Here we give the proof for the bounds on the trace 
and operator norm. The proof for the Frobenius norm is 
omitted, since it is similar and simpler. 

Let p S S^ . Them Theorem[51implies that there exists 
p & S, with PA = PA, such that: 



d-l 



1 



P- P 



N + d-V N + d 
Using the triangle inequality we have that 



\\P- P\\l < TT—i 7IIPII1 



^iPA^lB-p). (22) 



N + d-1' 

1 2(d-l) , , 

•\{pa®\b-p)\\i^ , \ , (23) 



N + d-1 



iV + d- 1' 



where in the last step we used once more the fact that 
Pa ® Is — P is separable (and, therefore, positive). Rela- 
tion (|16p is thus proven. 



5 



For the operator norm, let be the eigenvector 

corresponding to the maximum (minimum) eigenvalue of 
p — p. It follows that 

\\p-~p\\ca = max(tr{(/?-p)|M+)('u+|},tr{(p-p)|u_)(u_|}). 

(24) 

On the other hand, 



tr{(p - ~p)\u+){u+\} = j^^ir{p\u+){u+\] - 



N- 



i^tr{(pA ® Is - p)\u+){u+\} 



< 



d-l 
N+d- 



J, (25) 



and 



tr{(/5 - p)\u-){u-\} = -jJt^tr{p\u-){u^\} + 
+ NTd-MiPA®iB~p)\u-)(u^\}<^Jt^^. (26) 

The first part of the corollary has been proven. 

li pG and N > 2, then p can be seen to be PPT. 
Since the PPT criterion implies the reduction criterion 
[27L[28l|. we have that pa'^^b — p> 0. This observation, 
combined with the techniques used to derive the first set 
of relations, allows to prove the second one. □ 

The above corollaries can be reformulated as: 

Corollary 6. Suppose S{S) is a d-neighbor of the set of 
all separable states S in terms of the trace distance: 

Si6)''^'\J{aeS'\\\p-a\\<5} (27) 
pes 

(remember that is the set of all quantum states in 
Then, the following relations hold: 



SAT 



C S 




C S{gN) ^S[2[^ 



(28) 
(29) 



where the approximations are granted to hold in the limit 
iV>d> 1. 

This corollary suggests that the upper bounds for Sp 
converge quadratically faster than those for . In other 
words, if these bounds were optimum, then we would 
have proven that the additional PPT constrain gives the 
DPS criterion a quadratic speed-up. 

It is then natural to wonder if such bounds are indeed 
optimal. We will argue that at least the scaling of the 
upper bounds for is correct, i.e., fixing dA and d, 
the maximum possible entanglement robustness of any 
bipartite state pab arising from an N Bose symmetric 
extension scales with as 0{1/N). 

To see this, let N = 2K — 1, and consider the iV + 1 
bipartite state given by 



K K 



I* 



c 



\- J2 1 o^^>, 



(30) 



where Ck is a normalization factor. Define now pab = 
trBiv-i(|*^Biv)(*^Biv|). Clearly, pab G . Now, it 
can be shown that 



Pab 



K-1 
2{2K~ 1) 



(|00)(00| 



111X111 



K 



2{2K ~ 1) 



(|01> + |10»((01| + (10|). 



(31) 



The partially transposed operator has a negative 
eigenvalue —1/2{2K — 1) corresponding to the eigen- 
vector (|00) - |ll))/\/2, whose maximum Schmidt co- 
efficient is l/\/2- According to [12], this implies that 
R{pab) = l/i2K - 1) = 1/iV. The bound (O is, there- 
fore, tight for dA = d = 2. Since for any pair for Hilbert 
spaces Ha , 'Hb of dimensions greater than 1 we can em- 
bed the previous family of states in B{T-Ca ^Hb), it fol- 
lows that the optimal upper bound on the entanglement 
robustness of partial traces of Bose symmetric extensions 
must scale as 0{1/N). On the other hand, the bound 
(flS)) guarantees that the corresponding value for 5'^ at 
least scales as 0{1/N'^), i.e.. Theorem [3] allows to derive 
an upper bound for the entanglement robustness that 
decreases asymptotically faster than the optimal upper 
bound in the general Bose symmetric case. 

Note that the above considerations also allow us to 
obtain a dimension-dependent lower bound on the max- 
imum possible entanglement robustness R^,„ of a state 
in S'^. Following the lines of [38 



'sup 

^ consider the state 

(T = p^g , with PAB given by equation ([31]) . Clearly, 
a e 5^, with dA= dB = d = 2^. As -1/{2N) is the 
only negative eigenvalue of p^g and, therefore, the sum 
of its positive eigenvalues adds up to 1-1- 1/(2 A^), the 
negativity of a [39| (i.e., minus the sum of the negative 
eigenvalues of ct^^ ) can be seen equal to 



L(M-1)/2J 



M 
2j + l 



(1+2^) 



1 xAf-2i-l 



(27V)2J + 1 





^ 2Ar) 








2 


(1^ 


1 \M 
N > 


- 1 M 




2 


^ 2N' 



(32) 



where the last approximation is valid in the limit of 
large N. Since R{a) > Af{a) it follows that 

-R(cr) ^ 0{\og{d)/N). That is, for fixed dimension d, 
RZp satisfies 0(log(d)/7V) < Rg^ < 0{d/N). 
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IV. COMPUTATIONAL COMPLEXITY OF 

WMEM(5) 

In this section, we will analyze the consequences of the 
previous results on separability from the point of view of 
Computer Science. Actually, there are several different 
ways to describe the separability problem as a computa- 
tional problem Q . We chose to focus our attention in an 
approximated separability problem called the weak mem- 
bership problem of separability. This "promise" problem 
(as opposed to a "decision" problem) roughly consists on 
deciding the separability of a given state, but allowing an 
uncertainty parameterized by S. In this Section we will 
derive upper bounds on the time and space complexity 
when we attack this problem via the DPS criterion. 

The In-biased^ weak membership problem is defined 
as follows M : 



Definition 7. Weak membership problem of separability 
(WMEM(S)) 

Given a bipartite quantum state p £ and rational 6 > 
0, assert either that 



p G S{5) or 
p ^ S, 



(33) 
(34) 



where S{5) is a S neighbor of S, i.e., S{S) = {ct G 5^ 
a & S C S\\\a- a\\i<S}. 



In the above definition, ||cli||i = tr(\/ ojuj'^), the trace 
norm of the operator u, although, in principle, we could 
have chosen other norms or distance measures as an ac- 
curacy parameter. 

WMEM(S') is, thus, an approximation of the conven- 
tional separability problem in the sense that an algorithm 
solving WMEM(5) may assert equation ([55]) for a state 
Pab having just a small amount of entanglement. This 
approximated formalism is more practical than a non- 
approximated or exact formalism like EXACT-QSEP 6], 
because of the inevitable errors we incur in both numeri- 
cal and experimental studies, that should somehow be ac- 
counted for in our analysis of separability. A fair amount 
of effort has been devoted to the study of the time com- 
plexity of WMEM(5), the most remarkable result being 
that, if dA> ds, then WMEM(S') is NP-hard whenever 
1 /6 increases exponentially [40] or polynomially [4l| with 
respect to ds- 

We will now proceed to evaluate the time complex- 
ity of WMEM(5) when solved through the DPS crite- 
rion. First, following the discussion of Dohcrty et al. 
[lo| . can be characterized by a semidefinite program 

with ^(dim7i|yj^)^ — d\ free variables and a matrix 
of size {dim'Hgyjy^)dA on which we will impose the pos- 
itivity constraint. On the other hand, for Sp , the PPT 
constraint implies demanding positivity from an addi- 
tional matrix of size (dimTlgy^)^dA. Since the time- 
complexity of an SDP with m variables and of matrix 



size n is Oirn^n^) (with a small extra cost coming from 
an iteration of algorithms), the dominant factors for the 
asymptotic time-complexity of these tests can be written 
as 



Symmetric : d^(dimHgym)^ (35) 
PPT symmetric : d\{AiiTYn^^^f{dunnsym^fl?>&) 



where A^sym and -/Vppt are the sizes of the extensions 
needed to achieve a given accuracy parameter 5. 

Thus, at this stage, even though converges to S 
faster than 5^ , there still remains the possibility that the 
algorithm based on the sets {S^} is slower than the one 
based on the sets {-5''^}, because of the increase in time 
complexity that arises from imposing positivity on the 
partially transposed operator. The following calculation 
will rule out this possibility. 

From Eq. ([28]) of Corollary [6l we have that 



■Tr < {2-S)(dB-l) 

J^sym = g : 

\/2jd„-2,l 



(37) 



Taking into account that ~ d+0{d^/^) [2l[, the final 
expressions for upper bounds of the time complexity with 
respect to one method and the other are 



o\d\ 



o\d\ 



2e 
T 



Me 



for S 



N 



for S. 



N 



(38) 



where we just wrote the dominant (exponential) terms 
and omitted all polynomially growing terms. Note that 
the scaling law derived for the non PPT DPS criterion is 
valid as long as the optimal bounds on the trace distance 
to the set of separable states scale as ds/N . We conjec- 
ture that such is the case, although all our attempts to 
derive an analytical proof have failed so far. Under this 
assumption, the above formula thus shows that the crite- 
rion based on PPT BSEs indeed requires less steps than 
the one based on plain BSEs in order to solve WMEM(5^) 
for a given accuracy 5. 

The space complexity of both the plain DPS criterion 
and the PPT DPS criterion, though, is of the same type. 
This is because, although the PPT condition imposes (at 
least) a quadratic speedup in the speed of convergence, 
it also increases quadratically the size of the matrices 
involved in the SDP. Thus one effect cancels the other, 
and the size of the matrices needed in both cases to solve 
WMEM(5^) up to a given precision 5 is comparable for 
any value of ds. It follows that, according to our bounds, 
in some situations it may be more convenient not to use 
the PPT condition in order to save memory space. 
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Our experience with the DPS method suggests, how- 
ever, that this expectation is not reaUstic, but rather a 
consequence of the non optimahty of the bounds impUcit 
in Theorem [21 Actually, in practice, the algorithm based 
on PPT BSEs seems to have smaller space complexity 
than the one based on general BSEs. 

A big underestimation of the role of the PPT condition 
in the DPS criterion could also explain why the bound 
([38]) behaves much worse than the asymptotic expres- 
sions derived in Q for the performance of the 
algorithm conceived by loannou et al. for entanglement 
detection [ist- Indeed, as we will see, our bounds on 
the distance between the sets Sp and the set of separable 
states are far from optimal, at least for small values of 
cIa- Therefore, a more refined analysis could in principle 
end up with a different scaling law for this distance, that 
would eventually lead to a much better estimate of the 
time complexity of methods based in PPT BSEs. 



choose arbitrary states as a guess. 

Being a pure state, the efficiency of the protocol 
as a whole can be parametrized in terms of the average 
fidelity f: 

< / = ^p,tr(^'^Af,)tr(0,*,) < 1. (39) 

And the state estimation problem consists on determin- 
ing F, the maximum fidelity among all possible measure- 
and-prepare schemes {Mxt4>x)- Since F can be used 
as well to determine whether a given quantum channel 
can be simulated or not by an entanglement breaking 
channel, this problem is also referred to as the Quantum 
benchmark problem (Sll [s^ . [ssl . [s^ . [ssj . 

In [56| , it is explained how to map the SE problem into 
a linear optimization over the set S of separable states, 
via the relation 



V. APPROXIMATE ALGORITHMS FOR STATE 
ESTIMATION, MAXIMUM OUTPUT PURITY, 
AND GEOMETRIC MEASURE OF 
ENTANGLEMENT 

There are many relevant quantities in quantum infor- 
mation whose definition involves a linear optimization 
over a set of separable operators. The maximum average 
fidelity in state estimation problems [i^, [s^l , the output 
purity of a quantum channel [2^ or the geometric mea- 
sure of entanglement [l| are examples of such quantities. 
In order to compute these functions, we could think of 
an approximate algorithm that optimized over the sets 
or instead of 5, and it is easy to see that such 
an algorithm would give the correct answer in the limit 
of large N . 

So far, we have seen how Theorems and [3] can be 
used to derive bounds related to the separability prob- 
lem. In this Section we will show how to use these same 
theorems to bound the precision of the approximate lin- 
ear optimizations over the cone of separable operators 
mentioned above. 



A. State Estimation Problems 

In a general state estimation scenario, a source chooses 
with probability pi a virtual quantum state ^'^ that is 
encoded afterwards into another quantum state 4" J, to 
which we are given full access. The goal of the game is 
to measure our given state by means of a Positive Oper- 
ator Valued Measure (POVM) {Mx]x and thus obtain a 
classical value x that we will use to make a guess (j)x on 
the original state ^f^, which from now on we will assume 
to be pure. In conventional estimation theory, we usu- 
ally restrict the guess (j)x to be one of the original states 
[H, In this section, however, we will con- 

sider the more general setting in which we are allowed to 



F = max{tr(pAi3AAi3) : ^ab eS,AA^ I}, (40) 

where pab — J2iPi^i ® given by the particular 

SE problem. There it is also shown that any separable 
decomposition of the optimal operator Aab — J2x ^'^^ ® 
(j)x corresponds to the optimal strategy {Mx,<f>x)- 

Now, consider the sequence of optimization problems: 



F'' = max{tr(pABAAB) : Aab eS^,AA= I}, 

F^ = max{tr(p^BAAB) : Aab e , Aa = I}, (41) 

From (121), it is immediate that F^ > F'^ > F^ > ... > F, 
with limN^oc — F. An analogous property holds 
for the bounds F^ . Note that these maximizations are 
SDPs and therefore can be easily computed. 

Unfortunately, given limited computational (and spe- 
cially memory) resources, it is only possible to compute 
these bounds up to some index N. In spite of the asymp- 
totic convergence of the sequence, F^ or F^ could very 
well be far away from the actual solution of the problem. 
Is there any way to estimate the error of the truncation? 

Take Aab & S^{Sp) to be the operator that maxi- 
mizes equation (HIT) . Theorem ^ then implies that 
Aab, as defined by equation ^ ((UHl)), corresponds to 
a feasible state estimation strategy, since it is separable 
and Aa — Aa — I- Moreover, we can use the separa- 
ble decomposition of Aab that appears in Section IVII to 
express it as a measure- and-prepare protocol (Mx, 4>x)- 

The fidelities F^ or F^ associated to these strategies, 
although non trivial, will not be optimal in general, but 
they should provide a lower bound for F. From (|40p. it 
is easy to see that 



F 



N 



N 



N + d 



-F 



N 



1 



Fp^^ll- 



N + d' 

N 



2(d-lW P 



9n 



2(d- 1) 



(42) 
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7 8 9 10 11 



FIG. 1: Upper (squares) and lower (circles) bounds for the 
maximum fidelity _F as a function of A'^. The dashed line 
indicates the value of the exact solution, attained exactly by 
the PPT upper bounds on F from N = 2 and onwards. The 
minimum difference between the upper and lower bounds is 
of the order of 10~^ in both plots. 



Notice that both lower bounds asymptotically converge 
to F. That is, from the solutions of the semidefinite pro- 
grams (|4T1) it is possible to obtain a sequence of state esti- 
mation strategies that converges to the optimal measure- 
and-prepare scheme. 

To have a grasp on the efhciency of the method, con- 
sider the following state estimation problem: suppose we 
have a device that outputs two copies of one of the 4 
qubit states {|^'fc)}|.=i = {|0), |1), |-)} with equal 
probabilities. Our task is to estimate the state produced 
by the device. However, due to the environmental noise, 
once we are ready to measure the copies, those have de- 
graded into pk = (| |) = (1 - e) | *fc)(^'fc | + eI/2. 
The results for e = 0.3 are shown in Figure [1] for both 
the PPT and non PPT case and different values of N. 

We used the MATLAB package YALMIP ^ in com- 
bination with SeDuMi |la | to perform the numerical cal- 
culations. Note that the curve corresponding to the up- 
per bounds is constant, i.e., F^ — i^*^ — F* , for all 
M, N. This suggested that F* could be equal to F, 
the solution of the problem, although we did not ob- 



serve any rank loop in the matrices output by the com- 
puter. We thus had to force the rank loop to occur. 
Using rank minimization heuristics [i^ we checked for 
the existence of low rank PPT BSEs of Aab such that 
tT{AABPAB) > F* — S. Taking S = lO^'^, the computer 
returned a matrix with a rank loop, therefore proving the 
optimality of F* up to this precision. 

We performed a similar analysis for d = 3, this time 
considering the problem where a degraded copy of one of 
the states 



I*, 



cosif |0) + ,i„ f 



cos 



m , , , 



(43) 



(where i and j run from to 5) is sent to us with proba- 
bility 1/36 through a depolarizing channel p — > Q'^^-^''{p). 
In this case we were also able to force a rank loop in 
the PPT BSEs, so we again knew the optimal solution. 
Figure [5] illustrates our numerical results. 

Note that, in both cases, the lower bounds on the solu- 
tion behave very similarly as the upper bounds given by 
the DPS criterion, as long as we are considering the non 
PPT case. In the PPT case, however, our bounds prove 
to be terrible, since the second available upper bound ob- 
tained through the DPS criterion already seems to attain 
the optimal solution. We will discuss briefly this topic in 
Section HXl 

The main features of the practical performance of 
the DPS criterion have already been illustrated above. 
Therefore, in the following two problems we will just stick 
to analytical results. 



B. Maximal output purity of quantum channels 

Let Lu he a quantum channel. The maximal output 
purity 29] v of uj is defined as 



1^ = max||w(p)|joo, 
p 



(44) 



where the maximization is to be performed over all nor- 
malized quantum states p. 

At first sight this quantity may seem extremely non 
linear. We will show that, actually, (HH) can be reformu- 
lated as a linear optimization over the set of separable 
states. 

Denote by ^Iab the Choi operator corresponding to cj, 
i.e., Lu{p) = tr a{^ab • ® p)- It follows that 



1/ = max ||(jj(/9)||oo = maxtr{QAB • ct (g) p). 



(45) 



with a,p>0, tr(p) = tr(fT) 
Or, equivalently. 
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FIG. 2: Upper (squares) and lower (circles) bounds for the 
maximum fidelity _F as a function of A'^ in dimension 3. This 
time, the minimum difference between our lower bounds and 
the exact solution is around 0.03 (and it is attained in the 
non PPT case). 



C. Geometric entanglement of tripartite pure 
states 

Let \'^)abc be a pure tripartite state. A popular en- 
tanglement measure for this kind of systems is the so 
called geometric entanglement [s^, |3l[ (that some math- 
ematicians may recognize as the square of the e-norm 
[321), defined as 



max |(0^|(0s|((/.c||*Ai3c)r. (48) 

Notice, though, that, if we fix (j)A and (j)B, the state 0c 
maximizing the overlap will have to be proportional to 
{(PaH'PbW'^abc)- This overlap will be therefore equal to 



trc ( (<^A I ((As 1 1 ^' ABC ) ABC 1 1 0a) 1 ) 
= tr(pAB|0yi)(0A| «> |0b)(0b|), 



(49) 



where pab = trc(|^'ABc)(^ABc|)- It follows that E can 
also be reformulated as a linear optimization over 5, i.e.. 



E = max{tr(AABPAB) : A^b G S}. 



(50) 



As before, converging and decreasing sequences 
{E'^)n,{Ep)n of upper bounds on E can be derived 
via the DPS criterion, and Theorems [21 El allow us to ob- 
tain complementary increasing sequences of lower bounds 
{E^)N,iE^)N, given by 



E^ = 



~ 2{d-l) 



N jTiN I 1 \ 
P ' 2(d 



E^ + ^^\a. 



(51) 



Here Xa denotes the smallest eigenvalue of pA- 



VI. PROOF OF THEOREMS [H [1 



max{tr(ri^BAAB) : ^AB e S}. 



(46) 



As in the state estimation case, it is possible to define 
decreasing sequences {iy'^)N, {'^p)n of upper bounds on 
h' that converge asymptotically to the optimal output pu- 
rity of the channel. Using Theorems [21 and [3l together 
with the fact that tvBi^AB) = Ia, we have that there ex- 
ist sequences {C'p)n of lower bounds on given 
by 



The purpose of this section is to derive Theorems [2l [3l 
But first, a few words on notation. 

Given a unitary operator U, by \U) we will denote the 
state U\0). Also, for any permutation tt G Pat, K- G 
B{H^^) will represent the corresponding permutation 
operator. V alone must be understood as the SWAP 
operator acting over a bipartite system 7i^^, i.e., 



(52) 



.AT 



N 



.N 



N + d 



1 - 



dg 



1 



N + d' 



N 



9N 



2(d-i); p 2{d-i)' 



(47) 



To finish, Ti^m "^ill denote the symmetric subspace of 
H'^^ (the dimension ofH will be clear from the context). 

We will now proceed to proof Theorems [H [31 
The basic idea for both proofs is to notice that the orig- 
inal problem of finding a separable state [63 | Aab very 
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close to Aab from its BSE A^^^^ can be viewed as a 
probabilistic state estimation problem (ssj . 

Consider the following protocol, in which Alice plays a 
passive part: 



1. A copy of A^s' 
and Bob. 



is distributed to two parties, Alice 



2. Bob performs performs an incomplete measure- 
ment over , described by the POVM {M^ > 
0}a;, with Mx < I. As a result, he obtains either 
an outcome x or a FAIL message, indicating that 
his measurement has failed to produce an outcome. 

3. If Bob receives a FAIL message, then he makes 
it public. Otherwise, he prepares a state e 
BCHb), and both Alice and Bob would output 
the state ^B^i^^^AB^)^'''': ^^^^ probability Px = 

Px 

tr^jv (AfajA^jv). 

The state Alice and Bob will produce conditioned on 
a non FAIL message will be then given by 



Aab 



tTBN{MxA 



'EyPy 



(53) 



and is, therefore, a separable state. Moreover, since any 
entanglement breaking map can be decomposed as a mea- 
surement followed by the preparation of a state, this is 
the most general linear map we can apply over Ti.^^ in 
order to return a separable state A^^. 

But how to find a measure-and-prepare strategy for 
Bob such that Aab is close to Aab^ A possible scheme 
could be that Bob pretended that his subsystems are N 
identical copies of an unknown pure state, performed to- 
mography over each of these subsystems independently 
and then prepared a state consistent with the average 
values he would measure. This strategy should give good 
results in the particular case where A^^w can be approx- 
imated by a state of the form 



p{U)dUpu ® \U){U\ 



(54) 



However, supposing that the state had the form above, 
an even better strategy would be to allow Bob to perform 
collective measurements over his subsystems and then 
prepare the most convenient state. 

In conclusion. Bob should apply a POVM that allows 
him to efficiently identify the state C/|0)(0|C/^ out of N 
copies of it. Because in principle Bob has no a priori 
knowledge of p{U)dU, it is reasonable that he assumes 
that p{U)dU — dU, the Haar measure. 

In this particular case, the best state estimation strat- 
egy and the best probabilistic state estimation strategy 
coincide [s^. This implies that Bob should apply the 
POVM {\U){Uf'^dU}u and prepare the state \U){U\ 
whenever he gets the result U. Therefore, 



r _ / dUtiB^ (Ia <^ \U){U\^''+'Aab- ^ Ib) 

Jdmr{\u){u\^^ps«) ■ ^''^ 

To evaluate these integrals it is enough to notice that 
1. For any operator C, 

f dUU®^C{U^)®^ = c^V^, (56) 
for some coefficients c,r- In particular, 



dU\U){U\®'^ = 



{d~l)\N\Piy^ 
{N + d-l)\ 



{N + d-l)\ 



(57) 



2. Due to the fact that A^^w acts over Ha ® '^lym' 
for any tt G Pn+i, 



trB«{(A^Biv (g) Ib)Ia ® V^} = 

^(Aa® Ib, if T^{N + l)^N+l; 
I Aab^ otherwise. 

Finally, we arrive at the expression 



(58) 



Aab 



N 



N 



-A 



1 



AB 



N + d 



-.A. 



(59) 



We have just proven Theorem [2] 

The next step is to extend the previous ideas to account 
for the PPT condition, and a possible way is to modify 
the previous bipartite protocol to give Bob the ability to 
transpose part of his state before proceeding with any 
measure-and-prepare scheme. Suppose then that Bob 
partially transposes a partition _B', corresponding to half 
of Bob's systems in Aab" (we will take N even for sim- 
plicity). Following the previous arguments. Bob could 

pretend that he and Alice are sharing a state A^^' very 
similar to 



J p{U)dUpu {\U){U\ ® |C/*)(C/* 1)^^/^ (60) 

The benefits of this apparently useless step become evi- 
dent when we take into account the well established fact 
that it is easier to estimate a state from a copy and its 
complex conjugate than from two identical copies [33l.[3^. 
In the case of TV = 2, the optimal POVM has the form 
{U(g)U*\ip){ip\{U(S)U*ydU}, where \ip) is a linear combi- 
nation of |00) and 1^'"'") = J2i the (non normalized) 
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maximally entangled state. The optimal strategy for gen- 
eral N is not known, but we suggest the measurement 



(j)udU = {U ®U*f^/^m(j)\{U^ ® {U*yf^/^dU, (61) 

followed by the preparation of \U){U\. Here \(t>) is an 
arbitrary linear combination of the states [63j \(j)n} = 
|00)®«|^'+)^/2-", i.e., 



N/2 

E 

n=0 



Cn\4>n)- 



(62) 



Of course, applying the POVM (jju over is equiv- 

alent to apply the (non positive!) map associated to 
C/^^|0)((/)|^«'(?7^)®^/2 o^g^. o^j. s^ate A^siv. That way, 

we can use the same tricks employed in the computation 
of jnH). 

A fast way to perform these calculations is to notice 
that, for m > n, 



= |00>(00r"®(I®|0>(0|)^'"-"®V 



(63) 



Therefore, there exists a pair of permutations tt, tt' G P^r 
such that 



VM{^^\^^'Vl = |0)(0|«"+" ® (64) 
But Ia ® VjA^BN = Kab«^a ® A^B«, so 

trB«(A^B"IlA ® C/^^|0„)(0,nr-'(C/^)^^) - 

= trBiv(A^BivlA ® |f7)(f7|^'"+" (8)I<^^-"""). (65) 

In the end, we have that 



/ Ac\ Ac 

J^AB = 1 - d—r^_^ Aab + — --^Aa (E> Is, 
y c^BcJ c^Bc 

where A and B are square matrices given by 



(66) 



B 



(n+m)! 



(n+m)! 



"™ — („+m+d-l)! ' — (n+m+d)!' 

n,TO = 0,l,...,iV/2. (67) 

In case of odd N, we would make Bob partially trans- 
pose {N — l)/2 parts of his state and then use the fol- 
lowing (incomplete) POVM: 



(68) 



After the appropriate computations, we again arrive at 
expression ([3]), but the form of A and B changes to: 



B„ 



(n+m+iy. 7 
{n+m+dy. ' 



(n+m+1)! 
(ri+m+d-l-1)! ' 



n,m = 0, l,...,(iV-l)/2. 



(69) 



Obviously, in order to guarantee that Aab is close to 
Aab, it is in our interest to minimize the quantity 



c^Ac 
c^Bc 



(70) 



over all possible vectors c. Details on how to calculate 
the minimum of (j70[) , together with the expression of the 
optimal c can be found in Appendix [K\ The result is: 



min/Ar(c) 



2(d- 1) 



(71) 



This concludes the proof of Theorem [3l 

Notice that in both cases the given separable decom- 
position of the states Aab is continuous. Because of the 
presence of the Haar measure, however, via Design The- 
ory it is possible to arrive at an approximate [351] or exact 
[3g | finite separable decomposition for these operators. 



VII. EXTENSIONS TO MULTISEPARABILITY 

So far, we have only been considering separability in 
bipartite systems. In this section, we show that almost 
all the results we have derived can be easily extended 
to deal with separability in m-partite scenarios. More 
concretely, we will show how to generalize Theorems [2] 
and [3] to the multipartite case, since, as we have already 
seen, most of the other results are just corollaries of these 
two theorems. 

In this case, we will be interested in sets S'^ of states 
that derive from an N locally (PPT) Bose-symmetric 
extension 



37| 



Definition 8. TV locally Bose-symmetric extension 
Let Ai23... e B{Ti.i ® 7^2 ® ® ...) be a non negative 



operator. We will say that A]^2™3" 
A123... iff: 



.) is an N locally Bose symmetric extension of 



1. A12N31V 

2. tr2N -i-^N - 



> 0. 

i...(Ai2«3iv ) ~ A123.... 



3. A]^2™3™... independently Bose symmetric in sys- 
tems 2, 3, 4.... 

As before, in case such extension is PPT with respect 
to some partition, we will denote it as an N PPT locally 
Bose-symmetric extension. 

How close is A123... to the set of separable states? 
Consider a triseparable system, for instance, and sup- 
pose that we have an N locally Bose-symmetric extension 
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h.j^gN(jN for Kabc- In order to estimate the distance of 
Kabc to the set of triseparable states we could conceive a 
protocol where the state A^^jvj^n is distributed between 
Alice, Bob and Charlie. As before, Bob and Charlie could 
then independently apply probabilistic state estimation 
over their subsystems and prepare both a quantum state 
depending on their measurement outcomes. 

From what we already have, the derivation of the final 
expression of the triseparable state Kabc is straightfor- 
ward. Equation ([7]) describes the action of Bob's strat- 
egy over any bipartite state. Considering the partition 
AC^\B-^ , it follows that the resulting tripartite state 
after Bob performs state estimation will be: 



N 



-A 



1 



ABC 



N + di 



-A^c-jv (g) If 



(72) 



Now it is Charlie's turn. This time we will take the 
partition AB\C^ . The final result is that 



It is interesting, thus, to try to determine how good 
the PPT criterion is for entanglement detection alone, 
i.e., not in combination with Doherty et al.'s method. 
Here, through a very simple argument, we show what we 
believe is the first result in this direction after the seminal 
paper of the Horodeckis [58| . 

The main idea of our derivation stems from the fact 
that positivity under partial transposition is equivalent 
to separability in C'^ ® systems [l^- Suppose, then, 
that we have a PPT state pab S B{'Ha ^ T^-b), with 
> 3, and dB 1^ '2, and consider the (non normalized) 
state PAB given by 



PAB « / dUdWP^ ® PwPABPf, ® Pw^ (77) 



where dU and dW denote the Haar measures correspond- 
ing to S{dA) and SU{dB), respectively, and 



^ABC = (N+dB){N+dc) ^^BC + [N+dJ)(N+dc) ^^B ® Ic + 

-,hAC ®Ib + jwTd±N+d^^A ® Isc (73) 



' (N+dB)(N+dc) 

is a triseparable state. 

The generalization to more parties is immediate. In- 
voking again the definition of depolarizing channels ([3]), 
in m-partite separability the expression for A1234... would 
be 



P3 ^ uY,\k){k\U\P^ ^wY,\k){k\W^. (78) 



k=0 



fc=0 



It follows that PAB is a convex combination of un- 
normalized states pu,w 

EE P3 ^ P^PabP^ ® Pw^ with 
pu.w G -B(C^(8)C^). Notice, also, that each pu,w is PPT, 
since 



where 



A1234... = (Il(^r!(p.))(Ai234...), 

i=2 



N + d^ 



(74) 



(75) 



The corresponding expression for A123... when it arises 
from an N locally Bose-symmetric extension, PPT with 
respect to the partition 12 3 ...|2 SL^/^J 
is still dill), but this time 



{ph 



PwPabPu • 



p2 -^Tb 



— p-^ 



P^,PabPu<^Pw'>0- (79) 



Since PPT equals separability in C'^ (8) systems, it 
follows that each pu,w is separable, and so is pab, since 
by construction it is a convex combination of these states. 

It only rests to find an analytical expression for pab- 
Using the previous techniques it is straightforward to ar- 
rive at 

Theorem 9. Let pab e B{TiA®'HB) be a PPT normal- 
ized quantum state, with dA ^3,dB > 2. Then, for 



2(d. - 1) 



9n 



(76) 



VIII. THE POWER OF PPT ALONE 



The Peres-Horodecki criterion, aka the PPT (Positive 
Partial Transpose) criterion 2], is one of the most pop- 
ular existent criteria for entanglement detection. It is 
simple, it provides a very good approximation to the set 
of separable states in small dimensional cases and it usu- 
ally leads to analytical results when applied over families 
of quantum states. Actually, some entanglement mea- 
sures, like the negativity [3^ or the PPT entanglement 
robustness [J are based on the PPT condition. 



dA{dA-3) dB[dB~2) 

PA = -72 ^PB = —72 \ , (80) 



d\-l 



the state fl'-P^^ (g) H^P'^Hpab) is separable. 

Note that, in the particular case dA — S, dg — 2, pab = 

PAB- 

By simple application of the tools already developed, 
we end up with the following Corollary. 

Corollary 10. For any PPT state pab, with dA > 
^,dB > 2, 



Rg{pab) < j^idA + l){dB + !)-!, (81) 
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FIG. 3: Optimum bound on the global robustness of entan- 
glement R for generic states (dashed line), as opposed to the 
upper bound for PPT states (solid line). In this plot, we as- 
sume that dA = ds ~ d. Note that the new bound becomes 
trivial as soon as d > 9. 



and there exists a separable state a such that 



\PAB - CT||l 



< 2- 
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{dA + l){dB + l)' 



(82) 



To get an idea on how good these bounds are, have a 
look at Figure ini There the maximum possible global ro- 
bustness of entanglement of a x state is compared 
with our upper bound for PPT states. We see that, al- 
though our upper bound becomes useless for d > 9, it 
is very powerful in the small dimensional case. For in- 
stance, for C'^ X C'^ systems, the bound is equal to 1/3 as 
opposed to 2. This means that we would have to apply 
the non PPT version of the DPS method up to iV = 6 in 
order to characterize likewise the set of separable states. 



IX. CONCLUSION 

In this paper, we have studied the efficiency of the DPS 
criterion for entanglement detection. First, we showed 
that it is enough to subject the DPS states to some local 
noise in order to deprive them from their entanglement 
properties. It turned out that, while the minimal amount 
of noise necessary to turn an arbitrary state in into 
a separable state decreases as 0{\/N), the correspond- 
ing amount of noise needed to disentangle states in 
decreases at least as 0{\/N'^). We used these expres- 
sions to estimate the time complexity of both methods 
when applied to solve the Weak Membership Problem 
of Separability, and concluded that the PPT condition 
is worth imposing provided that the optimal bounds on 
the speed of convergence of the method based on plain 
BSEs scale as 0{d/N), as our own bounds suggest. We 
therefore hope to have shed some light on the question 



of how much the DPS criterion owes its strength to the 
PPT condition. 

We also derived bounds on the error we incur when we 
substitute the set of separable operators by or Sp 
in linear optimization problems, like the state estimation 
problem, the problem of determining the maximal out- 
put purity of an arbitrary quantum channel and the com- 
putation of the geometric entanglement. We performed 
numerical calculations of the first of these problems to 
test the accuracy of our analytical bounds. In order to 
compare our uncertainty with the actual solution of the 
problem, we developed a new technique that allows to 
prove in some cases the optimality of the DPS relax- 
ations. We observed that, although the bounds for the 
non PPT case seem to be very accurate, the bounds for 
the PPT case are too big when compared with reality. 

This disagreement between theory and practice may 
be explained in part by the fact that our bounds do not 
take into account the dimensionality of Alice's system, a 
crucial fact when dealing with the PPT constraint [5a |. 
For all we know, our PPT bounds could be exact in the 
limit cIa oo. Our intuition, nevertheless, is that better 
bounds could be found by applying linear maps over the 
initial state pab in order to obtain a separable state pab , 
as we did, but whose separable decomposition would be 
given by a non linear map, unlike in our examples. Ac- 
tually, we already used that approach in Section IVIIII to 
bound the entanglement of PPT states. That kind of 
schemes, together with state estimation considerations, 
may allow in the future to obtain such better bounds. 
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APPENDIX A: MINIMIZATION OF ([70]) 

Take N even. Then it can be checked that 



A — 



{d-iy. 



B™ = / dx. (Al) 

Combining this relation with (|70p . it follows that 



/(c) 



1 lolEnil^nX-l'i^-^n^-^r-'dx 



d~ 1 



lolEnllcnX-lH^-xr-'dx 



(A2) 
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That way, we can see the minimization of /(c) as a 
minimization over the set of all polynomials Qn/2{x) = 
^ c„a;" of degree N/2. Making the change of coordi- 
nates y — 2x — \ we find that the above minimization is 
equivalent to 



Q 



1 t-i\QN/2{yr{i-vY-'dy 

7n 2{d-l) |Q^/2(y)|2(l - yY-2ay 



(A3) 



where QN/2{y) is an arbitrary polynomial of order N/2. 
This problem can be solved by means of the Jacobi poly- 
nomials. 

The Jacobi polynomials P^^''^\y) are a complete set 
of functions orthogonal upon integration in the interval 
[-1, 1] under the weight {l + yY{l-y)°' [2l|. Now, define 
the normalized Jacobi polynomials Pniv) as 



Pniy) 



P^-'''\y) 



\P, 



(d-2,0)| 



(A4) 



with 



-2,0)11 ^ 




(d-2,0) 



{yWd-yY-^dy. (A5) 



It is clear that we can express any QN/2{y) as a linear 
combination of normalized Jacobi polynomials of order 
less or equal than N /2. That is. 



N/2 

QN/2{y) = X! ^^P"(y)' 



(A6) 



for some coefficients e„. Because of the orthogonality of 
the PnS, when we input this expression in the integral of 
the denominator, we end up with 

/' \QN/2iy)\'ii ~ yf-'dy ^ 

To calculate the integral on the numerator, we can 
make use of the recurrence relation 



Cm,n = an, if m = ft, 

Pn, if m = ?i + 1, 
7„, if TO = n - 1, 
elsewhere. 



(AlO) 



Now we will proceed to diagonalize C. 

Let A be an eigenvalue of C. This means that there 
exists a vector {wi}^''o such that 



(an - A)i;„ + PnVn+l + JnVn-1 = 0, 



(All) 



with WAr/2+1 = 0. 

Choose a real number yo and try the ansatz Vn — 
Pn{yo)- From (|A8|) . it is clear that Vn will satisfy (jAlip . 
provided that 



A = 1 - yo, 
PN/2+liyQ) = 0. 



(A12) 



That is, any root of the polynomial PN/2+i{y) corre- 
sponds to an eigenvalue of C . 



But PN/2+i{y) has N/2 -I- 1 simple roots [2l|, so all 
the eigenvalues of C are obtained using this strategy. It 
follows that 



min/Ar(c) 



1 



■ min{l - X : P^n+ii^) = 0} 



2(d-l) N/2+1 

(A13) 

Let us remark that this is not the first time the zeros of 
the Jacobi polynomials naturally appear in state estima- 
tion problems [STj . 

The expression for the case of odd N can be derived 
in an analogous way taking into account that, this time. 



• ""^^ ^1". dx. 



{d-i)\ 

n X{l-Xf- 



{d-2) 



-dx. (A14) 



(1 - y)Pn{y) = anPn{y) + PnPn + l{y) +lnPn~l{y), (A8) 

that holds for some coefficients an, Pmln, with 70 = 
and 7„+i = /3„ [21]. Invoking again the orthogonality of 
the Jacobi polynomials, we have that 

mill / (c) = min e^Ce, ( A9) 

c |el =1 2{d — 1) 

where C is an {N/2-^1) x (A^/2-|-l) tridiagonal hermitian 
matrix given by 



APPENDIX B: OPTIMALITY CRITERION 
(RANK LOOPS) 

For some problems involving linear optimizations over 
the set S, it may happen (see [s^) that a particular re- 
laxation of the problem F'^ turns out to coincide with 
F . In this appendix we will show how this optimality can 
sometimes be detected. 

We will take inspiration from optimality detection in 
other hierarchies of semidefinite programs that appear in 



15 



scientific literature. Consider the hierarchy of semidefi- 
nite programs used in |4^] for the calculation of the max- 
imal violation of linear Bell inequalities. There the op- 
timality of a relaxation is detected when the rank of the 
matrix generated by the computer is equal to that of 
some of its submatrices. Remarkably, we can find similar 
results in the hierarchies of semidefinite programs defined 
by Henrion and Lasserre to minimize real polynomials in 
a bounded region of M" ji^ . 

The corresponding result in this scenario is the follow- 
ing: 

Lemma 11. Let Aab" be a BSE of Aab, PPT with 
respect to the partition AB^'^\B^^^ . If 

rank{AAB'^) ^ max{ranfc(Ayi3R'), ranfc(A3iv-/f )} (Bl) 
then Aab a separable operator. 

Following we will say that A^i^jv presents a rank 
loop when it fulfills condition (jBl[) . 

The proof of Lemma [Tl] follows trivially from an old 
result by Horodecky et al. [46| : 

Theorem 12. Let pab be a PPT bipartite quantum 
state. If 

rank{pAB) < rank{pA), (B2) 
then pAB is a separable state. 

See [i^ for a proof. 

The possibility of finding a rank loop in practice in 
cases where the optimization over the set Sp coincides 
with the optimization over S should not be surprising. 
Note that any (finite dimensional) separable state 



can be expressed as a finite convex combination of prod- 
uct states, i.e., 

K 

Aab = '^PiPi ® \i>i){'^i\, with pi > 0,Vi, (B3) 

i=l 

with 7^ li'jXi^jl: for i ^ j. Now, consider the 

PPT Bose symmetric extension of Aab given by 

K 

1=1 

Clearly, as TV tends to infinity, the vectors 
become orthogonal. It follows that K* = 
limjv— >oo rank(Ayi5iv) exists and is equal to rank(pj). 
Being the rank a natural number, this implies that there 
is an M such that, for any N > M, rank(A^5Af) = 
rank(AA5]v) = K*. That is, for any finite dimensional 
separable state there exists a PPT Bose symmetric ex- 
tension with a rank loop. 

Of course, the fact that for any separable state pab 
there exists a PPT BSE with a rank loop does not mean 
that our computer is going to return such an extension. 
Note, though, that, if at the same time we set our com- 
puter to the task of finding PPT BSEs of pab we also 
demand a rank minimization of these matrices (i.e., we 
look for PPT BSEs with minimal rank), at some point 
we will find a rank loop. 

Unfortunately, rank minimization of positive semidefi- 
nite matrices with linear constraints is in general an NP- 
hard problem [3, [5^. There are, however, heuristics [g^l 
that have proven to be very efficient for solving small- 
scale problems (that is, for small d). 



M. B. Plenio and S. Virmani, Quant. Inf. Comp. 7, 1 [14] 
(2007). 

A. Peres, Phys. Rev. Lett. 77, 1413 (1996) [15] 
D. Bruss, J. Math. Phys. 43, 4237 (2002) 

B. Terhal J. Theor. Comp. Sci. 287, 313 (2002) [16] 
A. Sen, U. Sen, M. Lewenstein, and A. Sanpera, e-pnnt 
arXiv:quant-ph/050S032. [17] 
L.M. loannou. Quant. Inf. Comp., 7, 335 (2007) 

R. Horodecki, P. Horodecki, M. Horodecki, and K. [18] 
Horodecki, e-print arXiv:quant-ph/0702225. 

O. Giihne, G. Toth, e-prmt arXro;0811.2803. [19] 

A. C. Doherty, P. A. Parrilo, F. M. Spedalieri, Phys. Rev. 

Lett., 88, 187904 (2002). [20] 

A. C. Doherty, P. A. Parrilo, F. M. Spedalieri, Phys. Rev. 

A 69, 022308 (2004). [21] 

G. A. Raggio and R. F. Werner, Helvetica Physica Acta, 

62, 980 (1989). 

R. F. Werner, Lett. Math. Phys. 17, 359 (1989). [22] 
CM. Caves, C.A. Fuchs, and R. Schack, J. Math. Phys. [23] 
43, 4537 (2002). [24] 



L. Vandenberghe and S. Boyd, SIAM Review 38, 49 
(1996). 

J.F. Sturm, SeDuMi, a MATLAB toolbox for optimiza- 



tion over symmetric cones, http://sedumi.mcmaster.ca 
J. Eisert, P. Hyllus, O. Guehne, M. Curty, Phys. Rev. A 
70, 062317 (2004). 

F. G. S. L. Brandao, R. O. Vianna, Phys. Rev. Lett. 93, 
220503 (2004). 

F. Hulpke, D. Bruss, J. Phys. A: Math. Gen. 38, 5573 
(2005). 

R. Konig and R. Renner, J. Math. Phys. 46, 122102 
(2005) 

M. Christandl, R. Konig, G. Mitshison, and R. Renner, 
e-pnnt arXiv:quant-ph/0602130. 

M. Abramowitz and I. A. Stegun, Handbook of Math- 
ematical Functions with Formulas, Graphs and Mathe- 
matical Tables, New York: Dover Publications (1972). 

G. Vidal and R. Tarrach, Phys. Rev. A 59, 141 (1999). 
G. Vidal, R.F. Werner, Phys. Rev. A 65, 032314 (2002). 
M. Hayashi, D. Markham, M. Murao, M. Owari, S. Vir- 



16 



mani, Phys. Rev. Lett. 96, 040501 (2006). 
[25] D. Cavalcanti, Phys. Rev. A 73, 044302 (2006). 
[26] M. B. Plenio, F. G. S. L. Brandao, Nature, Physics 4, 

873-877 (2008). 

[27] M. Horodecki and P. Horodecki, Phys. Rev. A 59, 4206 
(1999). 

[28] N.J. Cerf, C. Adami and R.M. Gingrich, Phys. Rev. A 
60, 898 (1999). 

[29] G.G. Amosov, A.S. Holevo, and R.F. Werner, Probl. Inf. 
Transm. 36, 305 (2001). 

[30] A. Shimony, Ann. N. Y. Acad. Sci. 755, 675 (1995). 

[31] T. Wei and P. Goldbart, Phys. Rev. A. 68, 042307 (2003). 

[32] see e.g. A. Defant and K. Floret, Tensor Norms and Op- 
erator Ideals, North-Holland (1993). 

[33] J. Fiurasek, New J. Phys. 8, 192 (2006). 

[34] N. Gisin and S. Popescu, Phys. Rev. Lett. 83, 432 (1999). 

[35] A. Arnbainis and J. Emerson, Twenty-Second An- 
nual IEEE Conference on Computational Complexity 
(CCC07), pp. 129140 (2007). 

[36] A. Hayashi, T. Hashimoto, M. Horibe, Phys. Rev. A, 72, 
032325 (2006). 

[37] A. C. Doherty, P. A. Parrilo, F. M. Spedalieri, Phys. Rev. 
A, 71, 032333 (2005). 

[38] S. Beigi, P. W. Shor, e-pnnt arXi?;;0902.1806. 

[39] M.B. Plenio, Phys. Rev. Lett. 95, 090503 (2005), J. Eis- 
ert, PhD Thesis, Potsdam (2002). 

[40] L. Gurvits, Proceedings of the thirth fifth ACM sympo- 
sium on Theory of computing, 10, New York, ACM Press 
(2003). 

[41] Sevag Gharibian, e-print arA'i?;;0810.4507. 

[42] L. M. loannou, B. C. Travaglione, D. C. Cheung, and A. 

K. Ekert, Phys. Rev. A, 70, 060303(R) (2004). 
[43] L. M. loannou and B. C. Travaglione, Phys. Rev. A, 73, 

052314 (2006). 

[44] M. Navascues, S. Pironio and A. Acin, New J. Phys. 10, 

7, 073013 (2008). 
[45] D. Henrion and J. B. Lasserre, IEEE Trans. Aut. Contr. 

51, 192 (2006). 

[46] P. Horodecki, M. Lewenstein, G. Vidal and I. Cirac, Phys. 

Rev. A 62, 032310 (2000). 
[47] J. Lofberg, YALMIP : A Toolbox for Modeling and 

Optimization in MATLAB, http: / /control. ee.ethz.chT^ 

joloef/yalmip.php. 



[48] R. Orsi, LMIRank: software for rank constrained 

LMI problems, http://users.rsise. anu .edu.au/' 

robert /Imirank / . 
[49] C.W. Helstrom, Quantum Detection and Estimation 

Theory, Academic Press, New York, (1976). 
[50] A.S. Holevo, Probabilistic and Statistical Aspects of 

Quantum Theory, North-Holland, Amsterdam, (1982). 
[51] K. Hammerer, M.M. Wolf, E.S. Polzik and J.I. Cirac, 

Phys. Rev. Lett. 94, 150503 (2005) 
[52] A. Serafini, O.C.O. Dahlsten and M.B. Plenio, Phys. Rev. 

Lett. 98, 170501 (2007) 
[53] G. Adesso and G. Chiribella, Phys. Rev. Lett. 100, 

170503 (2008) 

[54] J. Calsamiglia, M. Aspachs, R. Munoz-Tapia and E. 

Bagan, e-print arXiv:0807.5126. 
[55] M. Owari, M.B. Plenio, E.S. Polzik, A. Serafini and M.M. 

Wolf, New J. Phys., 10, 113014 (2008). 
[56] M. Navascues, Phys. Rev. Lett. 100, 070503 (2008). 
[57] P. Rapcan, J. Calsamiglia, R. Munoz-Tapia, E. Bagan, 

V. Buzek, e-print arXiT;;0708.1086. 
[58] M. Horodecki, P. Horodecki and R. Horodecki, Physics 

Letters A, 223,1 (1996). 
[59] J. David, Algorimth analysis for robust controllers, PhD 

thesis, Kat. Univ. Leuven, ESAT, 3001, Leuven, Belgium 

(1994). 

[60] M. Fazel, PhD Thesis, Stanford University (2002). 

[61] To understand why, write p as a convex combination of 
product states, i.e., p = '^Pi\ui){ui\®\vi){vi\. Then, pA® 
I - p = TjVi\ui){ui \ (gi (I - \vi){vi\). That is, pA®l - p ^ 
Pa (g I — p is a separable operator. 

[62] Along this proof, we will consider the operator Aas to 
be a quantum state rather than a quantum operator, i.e., 
to be normalized. It is easy to see that, if ((7} and (|10p 
are separable states when Aas is normalized, the very 
same expressions have to lead to separable operators with 
tr(AAB) = tr(AA_g) if Aas is not normalized. 

[63] Note that a further symmetrization of these states over 
the particles in B' and B 

B' would be more intuitive, but irrelevant, since the sup- 
port of the state A^'„ is in Ha ® H^^m ® ^^m- That 
is, such symmetrization is automatically performed when 
we apply this POVM over A^^'„ . 



